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Abstract 


Land-atmosphere interactions play a critical role in determining the diurnal evolution of both 
planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of 
coupling between the land surface and PBL in numerical weather prediction and climate models remains 
largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a 
range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of 
Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process 
understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a 
framework for diagnosing local land-atmosphere coupling is presented using a coupled mesoscale model 
with a suite of PBL and land surface model (LSM) options along with observations during field 
experiments in the U. S. Southern Great Plains. Specifically, the Weather Research and Forecasting 
(WRF) model has been coupled to the Land Information System (LIS), which provides a flexible and 
high-resolution representation and initialization of land surface physics and states. Within this 
framework, the coupling established by each pairing of the available PBL schemes in WRF with the 
LSMs in LIS is evaluated in terms of the diurnal temperature and humidity evolution in the mixed layer. 
The co-evolution of these variables and the convective PBL is sensitive to and, in fact, integrative of the 
dominant processes that govern the PBL budget, which are synthesized through the use of mixing 
diagrams. Results show how the sensitivity of land-atmosphere interactions to the specific choice of 
PBL scheme and LSM varies across surface moisture regimes and can be quantified and evaluated 
against observations. As such, this methodology provides a potential pathway to study factors 
controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which will serve as a 
testbed for future experiments to evaluate coupling diagnostics within the community. 
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1. Introduction 


Land-atmosphere (L-A) interactions and coupling remain weak links in current 
observational and modeling approaches to understanding and predicting the Earth- Atmosphere 
system. The degree to which the land impacts the atmosphere (and vice-versa) is difficult to 
quantify given the disparate resolutions and complexities of land surface and atmospheric models 
and lack of comprehensive observations at the process level (Betts et al. 1996; Angevine 1999; 
Entekhabi et al. 1999; Betts 2000; Cheng and Steenburgh 2005; Gu et al. 2006). However, the 
convective planetary boundary layer (PBL) serves as a short-term memory of land surface 
processes (through the integration of regional surface fluxes on diurnal scales), and is therefore 
diagnostic of the surface energy balance. Further, the equilibrium established between land 
surface and mixed layer fluxes and states in a growing PBL is a function of the degree of 
coupling and the impact of feedbacks within the L-A system (Pan and Mahrt 1987; Oke 1987; 
Stull 1988; Diak 1990; Garratt 1992; Dolman et al. 1997; Peters-Lidard and Davis 2000; Cleugh 
et al. 2004; Betts and Viterbo 2005). As such, knowledge of temperature and moisture evolution 
in the PBL can be instrumental in estimating surface fluxes and properties across regional scales 
as well as quantifying and improving L-A representations in coupled models. 

Recent efforts to better understand and quantify the nature and processes involved in L-A 
interactions have focused on advancing the theory and formulation of their complex behavior 
and feedbacks (Sorbjan 1995; Steeneveld et al. 2006), deriving and exploiting relationships 
among observed L-A properties (Eltahir 1998; Santanello et al. 2005, 2007), and assessing these 
interactions using coupled models (Margulis and Entekhabi 2001; Barros and Hwu 2002; Ek and 
Holtslag 2004; Santanello et al. 2007). While progress has been made in identifying individual 
processes and feedback loops for a particular location or model, a comprehensive approach to 
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diagnosing the full nature of L-A coupling that can be applied to models and evaluated against 
observations has yet to be developed. The need for such a framework will only become more 
critical to ensure that advances in measurement technologies such as satellite remote sensing of 
the land surface and PBL are properly incorporated into L-A studies and models (e.g. data 
assimilation). 

A relatively simple but untested approach by Betts (1992) to quantifying heat and 
moisture budgets in the PBL is based upon a vector representation of the diurnal evolution of 
temperature and humidity. Application of this ‘mixing diagram’ theory to models and 
observations would offer the ability to perform a robust evaluation of L-A interactions with 
minimal inputs due to the integrative nature of the mixed layer on diurnal time scales. Ideally, 
this approach should be tested using a coupled, high-resolution, mesoscale model with flexible 
land surface and PBL schemes, thereby allowing the variation in L-A coupling among different 
formulations versus that observed to be evaluated. 

With these issues in mind, this paper defines a methodology to quantify local L-A 
coupling and the various components and feedbacks therein. Section 2 presents an overview of 
recent progress in L-A research and the complexities of the governing processes and feedbacks, 
including the growing need for studies of L-A coupling at the local (regional) scale, and 
describes the mixing diagram approach that is adopted and extended in this study. The coupled 
regional model, land surface models, and PBL schemes used in the experiments are highlighted 
in Section 3 along with detailed infonnation on the sites, case studies, and associated 
observations. Results and analyses of the mixing diagram approach applied to these experiments 
are presented in Section 4, followed by a discussion of the greater applicability and limitations of 
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this methodology in Section 5. Finally, Section 6 presents the conclusions and summary of 
current and future work related to L-A research. 

2. Background 

a. Motivation for Studying L-A Coupling 

The need for improved understanding, estimation, and prediction of L-A interactions and 
feedbacks has been growing significantly over the last decade (Jacobs and DeBruin 1992; 
Entekhabi and Brubaker 1995; Kim and Entekhabi 1998; Entekhabi et al. 1999; Liu et al. 2003, 
2004, 2005; Medeiros et al. 2005; Dirmeyer et al. 2004). During this time, offline (uncoupled) 
land surface models (LSMs) have grown in complexity and diversity, while the applicability of 
offline model evaluations, such as those performed during the Project for Intercomparison of 
Land-surface Parameterization Schemes (PILPS) experiments (Henderson-Sellers et al. 1993), 
may be severely limited to the omission of L-A interactions and feedbacks (e.g. Liu et al. 2003, 
2004, 2005). At the same time, LSMs coupled to atmospheric models are often highly tuned to 
each other without regard for the degree and accuracy of coupling between the L-A schemes or 
the impact of feedbacks. In both instances, our ability to diagnose and quantify these 
interactions is lacking, and needs to be improved by evaluating the best available PBL and land 
surface data in the context of the mixed layer evolution and equilibrium established through their 
interactions and feedbacks, and how this compares to what is simulated in our models. 

The Global Land/ Atmosphere System Study (GLASS), part of the Global Energy and 
Water Cycle Experiment (GEWEX; www.gewex.org) program, was designed to serve as an 
interface between the land-surface community and efforts to observe, understand, and model the 
hydrological cycle and energy fluxes in the Earth-Atmosphere system. GLASS is composed of 
four actions that support the intercomparison and advancement of a) large scale offline, b) large 
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scale coupled, c) local scale offline, and d) local scale coupled models. Currently, each of these 
is being addressed through organized, community-wide modeling studies [a) Global Soil 
Wetness Project (GSWP-2; Dirmeyer et al. 1999), b) Global Land Atmosphere Coupling 
Experiment (GLACE; Koster et al. 2002), and c) PILPS], with the exception of local scale 
coupled studies (hereafter referred to as ‘LoCo’). There have been community workshops and 
tentative implementation plans (e.g. van den Hurk et al. 2005) for a LoCo initiative, but to date 
an organized effort and lead for LoCo is still under development. 

Recent GLASS-related studies have highlighted the importance and difficulties in 
understanding the complexity and impacts of L-A interactions. For example, the GLACE study, 
while looking from a global perspective on non-local impacts of the land on the atmosphere and 
vice-versa, show that there are regions of highly coupled L-A enviromnents that impact 
precipitation patterns and cycling (Koster et al. 2004; Lawrence and Slingo 2005). These 
hotspots are therefore likely to be of interest to LoCo. Similarly, the GEWEX Atmospheric 
Boundary Layer Study (GABLS) community has attempted to isolate and intercompare an array 
of PBL schemes (i.e. single-column models) while controlling atmospheric and land surface 
boundary forcing, but in the process has shown the importance of accounting for a variable land 
surface that is fully interactive with the PBL (Holtslag et al. 2007). As a result, the GLACE and 
GABLS communities also have a vested interest in a LoCo action. 

From outside the GEWEX community, there have been a host of studies focused on a 
variety of individual L-A processes and feedbacks that call for further study of LoCo in a 
comprehensive and quantitative manner (e.g. Brubaker and Entekhabi 1996; Kim and Entekhabi 
1998; Berbery et al. 2003; Findell and Eltahir 2003ab). For example, Cheng and Steenburgh 
(2005) and Gu et al. (2006) demonstrate that the large variability of coupling in models and the 
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lack of quantitative understanding of the relevant processes need to be addressed including both 
direct and indirect (i.e. feedback) effects. Betts (2000) and Betts and Viterbo (2005) show that 
the L-A coupling is also critical in global model reanalysis data, but also stress that the critical 
processes and relationships that determine model evolution and equilibrium lie on the local scale 
and lack sufficient understanding and representation in models of all scales. 

Inherent in the ability to accurately simulate L-A interactions in coupled models is the 
engineering of the coupling itself in tenns of model design and variable passing from land 
surface to surface layer and PBL schemes, and vice-versa. Polcher et al. (1998) and Best et al. 
(2004) have proposed a generalized coupling design in this regard, but highlight the complexities 
involved in time-stepping of the coupled variables, and in specifying the blending height (i.e. 
tiling; Molod et al. 2004) and surface layer (Chen et al. 1997a) within each model. Overall, 
coupling design remains a largely model-dependent decision based on ease of implementation, 
rather than that which provides the most accurate representation of the L-A processes. 

Because we do not yet know the full nature of L-A coupling in models and acknowledge 
that the physics within LSMs is incomplete in many respects, there have been numerous efforts 
to calibrate model parameters to improve simulations. For example, Liu et al. (2003, 2004, and 
2005) demonstrated the ability to optimize large sets of L-A parameters, but that there are large 
differences in the offline and coupled cases. Unfortunately, as shown by Hogue et al. (2005), 
such extensive parameter calibrations are also completely model and site-dependent, and as a 
result do not tell us anything about the true nature of L-A coupling, its quantification, or how 
accurately it is represented in each model. 

Another important motivation for further understanding L-A coupling lies in its direct 
impact on data assimilation of PBL and land surface states. There have been numerous efforts to 
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assimilate screen-level observations in offline (Rhodin et al. 1999; Hess 2001), single-column 
(Hacker and Snyder 2005; Hacker and Rostkier-Edelstein 2007), and fully coupled (Seuffert et 
al. 2004) models. However, a great deal of testing of land data assimilation of soil moisture, 
surface temperature, and snow has been performed for offline models (Rodell and Houser 2004; 
Reichle et al. 2007; Bosilovich et al. 2007), which lack L-A interactions and feedbacks that 
otherwise would impact the assimilation results in coupled mode. As efforts to assimilate new 
remote sensing data increase along with the complexity of assimilation techniques, the manner in 
which the land and atmosphere are coupled as well as the strength of feedbacks becomes critical 
to the process. 

Recent attempts to further understand and quantify L-A interactions have identified key 
properties, relationships, and feedback mechanisms using targeted modeling and observational 
approaches. Ek and Holtslag (2004) have derived a formulation for relative humidity tendency 
at the top of the PBL that aids in identifying moistening and drying regimes, and incorporates the 
full set of L-A processes (and inherent feedbacks) governing PBL evolution. Similarly, Barros 
and Hwu (2002) investigated the L-A interactions that determine summer rainfall patterns using 
a mesoscale model. While their approach is slightly broader than that of Ek and Holtslag in 
spatial (synoptic) and temporal (daily-weekly) scales, they identify different timescales of 
feedbacks, as evidenced in the relationship of near-surface relative humidity to surface Bowen 
ratio, which represent moist and dry dominated processes. From a more local perspective, 
Santanello et al. (2005, 2007) demonstrate the utility of readily observable properties of the PBL 
and land surface (e.g. PBL height, diurnal 2m-temperature change, soil moisture, and 
atmospheric stability) and the strong relationships therein. Their results show that these 
properties are integrative of L-A interactions and feedbacks, and have since been supported by 
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studies using a coupled regional model (Desai et al. 2006; Reen et al. 2006) that highlight the 
importance of LoCo in atmospheric modeling. 
b. Mixing Diagram Approach 

1) METHODOLOGY AND APPLICABILITY 

The previous section details the recent progress, difficulties, and need for further research 
to quantify LoCo in observations and models. It is apparent from these studies that in order for a 
robust methodology to diagnose coupling to be effective and useful to the community, it must be 
comprehensive and integrative of L-A processes and feedbacks; while at the same time able to be 
implemented using easily observed and understood properties of the system. 

An approach that may satisfy these requirements for diurnal time scales is the concept of 
vector representation of heat and moisture (energy) budgets, as introduced by Betts (1984, 1992) 
in the fonn of ‘mixing diagrams’. This conservative variable approach relates the diurnal 
evolution of 2m-specific humidity (q) and potential temperature (0) to the surface and mixed 
layer energy balance and, in effect, the diurnal equilibrium established by L-A interactions. The 
daytime variability of 6 and q is sensitive to and integrative of the dominant processes involved 
in LoCo, and when plotted in energy space ( Lq vs. C P Q) can be used to quantify these processes. 

Figure 1 presents a mixing diagram of the temporal change in Lq vs. C p 6 as generated by 
a mesoscale model and representative of conditions during June 2002 at a point in Oklahoma. 

For a full derivation and discussion of this theory we refer the reader to Betts (1992). The 
temporal change from U to fy is fully described by vector components that represent the fluxes of 
heat and moisture from the land surface (V S f C ) and the top of the PBL (i.e. entrainment; V en t). 
These two vectors have a slope exactly equal to the Bowen ratio of the surface ( /fy ) and 
entrainment (f3 e nt), respectively. Their magnitude of their components, in terms of C P A9 and 
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LAq, are proportional to the fluxes of heat (//) and moisture (LE) of each, respectively. For 
example, the magnitude of the surface vector component in the y-direction (heat) is as follows, 


C p AG 


sfc 


H sfc M 

p m ■ PBLH 


( 1 ) 


which can be calculated from the mean surface sensible heat flux (I Lf) over the timestep (At = tf 
- 1,), the specific heat (C p ), the density in the mixed layer (p,„), and the mean height of the mixed 
layer ( PBLH) over At. PBLH therefore represents the volume by which the fluxes are scaled and 
can be defined as a suitable average of the depth of the PBL over the timestep (see Section 2b2) 
below). Note also that in Eq. 1 we have assumed that 6 and T are nearly equivalent at 2 meters 
(i.e. (0/T)-H s f c = H sfc ). 

Once the surface vector (V S f C ) is computed, the heat entrainment flux (H ent ) can be 
calculated from the component (C p AO enl ) of the residual vector that connects V S f C to the final 
values of C p 0 and Lq at tf. The analogous formulation for components in the x-direction ( LAq s f c , 
LAq ent ) is then used to derive to the fluxes of moisture at the surface (evaporation; LE s f c ) and top 
of the PBL (dry air entrainment; LE ent ). Figure 1 also shows the hourly time series of 0 and q, 
H ent , and PBLH in a more traditional sense, and in this example highlights the mid-late afternoon 
decrease in moisture and rise of the PBL associated with increased dry air entrainment. 

Having derived the slope (/f,,,) and magnitude ( C p A9 ent ) of V en t, the entrainment ratio (A R 
= HeJHsfc), defined as the proportion of heat input to the PBL from entrainment to that of the 
surface, can be easily quantified. Typically, estimates of A R are difficult to acquire and there has 
been little consensus as to what the value should be based on different empirical studies (Kustas 
and Brutsaert 1987; Betts and Ball 1994). The mixing diagram approach is extremely valuable 
in this regard, considering that difficulties in observing and measuring fluxes at the top of the 
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PBL greatly limit efforts to close energy budgets in the PBL (e.g., Peters-Lidard and Davis 2000; 
Santanello et al. 2005). Similarly, the entrainment of moisture (typically negative due to drier air 
in the free atmosphere) is easily quantified using this methodology. Therefore we find it useful 
to define separately a heat and moisture entrainment ratio (A/, and A/ e , respectively) when 
discussing the components of the energy budgets derived from mixing diagrams. 

Betts (1992), Betts and Ball (1994) and Betts et al. (1996) provide the foundation for this 
approach and apply it to empirical data from short-term field experiments. They show how the 
diurnal evolution of q and 0 strongly reflects conditions and processes at the land surface (soil 
moisture, evaporation) and the top of the PBL (entrainment) as the theory suggests. A 
qualitative example is provided by Betts et al. (1996; their Fig. 8) that compare mixing diagrams 
for three days exhibiting different soil moisture conditions in Kansas. Although the vector and 
flux components are not explicitly calculated, a visual examination of these curves indicates the 
different evolution of temperature and humidity for dry, intennediate, and wet soils. In 
particular, the impact of entrainment is most notably visible in the mid-late day drying of the 
mixed-layer, and is most obvious for dry surface conditions. In contrast, wet surfaces evaporate 
more freely throughout the day, moistening the shallower PBL and reducing the magnitude of 
entrainment, and also evidenced in the diagrams. 

2) RELATION OF STEP-WISE INTEGRAL TO DIURNAL APPROXIMATION 

While the mixing diagram approach is relatively straightforward, it is important to 
address its applicability to the diurnal cycle as proposed here (Fig. 1). As indicated by Eq. 1, the 
derived fluxes and metrics are rather sensitive to the values of H s f c , PBLH, and At, particularly 
early and late in the day when surface fluxes are small and the PBL is transitioning between 
nocturnal and convective behavior. To address this issue, we compared the results of applying 
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the mixing diagram approach hourly (or step-wise, with values of 0, q, H s f c , and PBLH varying at 
each hourly timestep) with those from a daytime mean (using the initial and final values of 9 and 
q, and the mean values of H s / C and PBLH between t, and tj) calculation of the fluxes and ratios 
derived from the vector components. This was performed using output from the coupled 
mesoscale model employed in this study at the sites presented in Section 4a. 

Despite lacking the finer temporal variability of the step-wise integration, the daytime 
mean fluxes (R 2 > 0.98) and entrainment ratios (R 2 > 0.97) at each site were found to correlate 
strongly with those calculated using hourly values averaged across the time period. While the 
step-wise approach is assumed to be more accurate, the small cost (MSE ~ 25.0 Wm’ 2 ) of simply 
using only the initial and final points with the mean flux and PBLH throughout the day suggests 
the derived fluxes scale nearly linearly. Inspection of the variability of hourly flux and 
entrainment ratios throughout the day indicate that the highly sensitive time periods actually 
balance out, with higher entrainment ratios in the morning when surface fluxes are low (and 
therefore entrainment large in proportion) and lower ratios in the late afternoon when PBLH is 
large relative to surface fluxes. This results in the average fluxes and ratios over the entire 
period to be adequately represented by the daytime mean approach. 

By using the theory outlined above, the exact nature of the L-A processes controlling 
PBL evolution, heat and moisture budgets, and their critical feedbacks can be evaluated in 
coupled models and compared against observations. It should be noted that the temporal 
dynamics (e.g. hourly evolution) of the fluxes and metrics can also be instructive in 
understanding the relative importance of PBL versus land surface processes in generating the 
daily mean L-A equilibrium. However, in terms of acquiring bulk infonnation on L-A coupling 
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throughout the day and being able to evaluate it using observations, the daytime mean approach 
remains more sensible and valuable. 

3. Model and Site Description 

a. WRF and LIS-WRF 

The Advanced Research version of the Weather Research and Forecasting (WRF-ARW) 
model (Michalakes et al. 2001) is a state of the art mesoscale numerical weather prediction 
system. Derived from the Fifth-Generation NCAR/Penn State Mesoscale Model (MM5; Anthes 
and Warner 1978), WRF-ARW has been designated as the community model for atmospheric 
research and operational prediction and is ideal for regional simulations on the order of 1-7 days. 
The model has a Eulerian mass dynamical core and includes a wide array of radiation, 
microphysics, and PBL options as well as 2-way nesting and variational data assimilation 
capabilities. 

To serve as a testbed for LoCo diagnostics, WRF-ARW (Versions 2. 1.2/2. 2) has been 
coupled to NASA’s Land Information System (LIS; Versions 4. 2/5.0) by Kumar et al. (2008). 
LIS consists of a suite of LSMs and provides a flexible and high-resolution representation of 
land surface physics and states which are directly coupled to the atmosphere (hereafter referred 
to as LIS-WRF). The advantages of coupling LIS to WRF-ARW include the ability to spin-up 
land surface conditions on a common grid from which to initialize the regional model, flexible 
and high-resolution soil and vegetation representation, additional choices of LSMs that will 
continue to grow, and various plug-in options such as land data assimilation and parameter 
estimation. LIS-WRF has been tested extensively thus far over the U.S. Southern Great Plains 
(SGP), Florida, the Gulf of Mexico, and Korea. 
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The LSMs employed in LIS for this study are the Noah LSM (Noah; Ek et al. 2003), and 
the Community Land Model Version 2.0 (CLM; Dai et al. 2003). Both models dynamically 
predict water and energy fluxes and states at the land surface, but vary in specific 
parameterizations and representation of soil and vegetation properties and physics. For example, 
the Noah LSM solves moisture and heat transport through 4 discrete soil layers while CLM 
solves for 10 layers. In addition, treatment of vegetation properties (such as Leaf Area Index and 
Vegetation Fraction) and canopy fluxes differ between the two FSMs. The Noah model 
employed in this study is Version 2.7.1 and is identical to the version of Noah packaged in the 
original version of WRF-ARW. CFM is unique to LIS-WRF, and also serves as the land model 
for NCAR’s coupled Community Climate System Model (CCSM). As such, these two LSMs are 
widely-used and capture the range in complexity (layering and vegetation physics) and 
application (mesoscale to global climate model) of schemes evaluated during the PILPS 
experiments. 

There are three options for PBL schemes in WRF-ARW, all of which are employed in 
this study using LIS-WRF. The Medium-Range Forecast (MRF; Hong and Pan 1996) scheme is 
based on non- local-K theory (Troen and Mahrt 1986) mixing in the convective PBL, where the 
diffusion and depth of the PBL are a function of the Richardson number (Ri cr ). The Yonsei 
University (YSU; Hong et al. 2006) scheme, based on the MRF, is also a non-local K theory 
implementation, but includes explicit treatment of entrainment and counter gradient fluxes. 
Finally, the Mellor-Yamada-Janjic (MYJ; Janjic 2001) scheme is the most complex PBF scheme, 
and employs nonsingular level 2.5 turbulent kinetic energy (TKE) closure (from Mellor and 
Yamada 1982) with local-K vertical mixing. In the MYJ scheme, the length scale is a function 
of TKE, buoyancy, and shear, and the PBF Height is diagnosed based on TKE production. 
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Overall, these three PBL schemes span the range in complexity (1 st order to TKE) and 
application (single column to full 3-D) of those employed in the GABLS experiments. We 
therefore expect the results to encapsulate the wide range of coupling possible between LSMs 
and PBLs participating in PILPS and GABLS. 

To address LoCo under the LIS-WRF framework, simulations were perfonned using the 
Noah and CLM LSMs with the MRF, YSU, and MYJ PBLs, for a total of 6 different 
combinations of L-A coupling (the remainder of the LIS-WRF setup is identical for each). The 
results of each experiment is then evaluated using the mixing diagram approach described above, 
where the processes and feedbacks generated by each LSM-PBL pair are quantified over the 
course of the day for different locations and conditions and compared with observations. 
b. LIS-WRF Experimental Design 

As shown by Koster et al. (2004) and others, the SGP region is a hotspot for L-A 
coupling in tenns of the strength of interactions and feedbacks. Because of this, and the wealth 
and record of observational data from the Atmospheric and Radiation Measurement (ARM) 
testbed located in the region (ARM-SGP), numerous intensive field campaigns have been 
conducted in the region that have augmented the instrument and data quality even further. For 
this study, simulations have been performed for the International FLO Project in June 2002 
(IHOP-02; Weckworth et al. 2004). During IHOP-02, we chose to focus on 36-hour simulations 
beginning at 00Z on June 6 and June 12, as they represent a clear-sky ‘golden’ day, and an 
unstable day with spatially heterogeneous convection, respectively. 

LIS-WRF simulations were perfonned over a large single domain (900x600), centered 
over the Oklahoma and Kansas border near the ARM-SGP central facility (CF) at Lamont, OK, 
with a horizontal resolution of 1 km and timestep of 5 s. The vertical grid of WRF-ARW uses an 
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r| level-formulation, and was specified with 45 levels and a lowest model level of 42 meters. 
Notably, there are approximately 21 model levels below 600 mb (4-5 km; the maximum depth of 
the mixed layer) which enables PBL structure and processes to be highly-resolved. Other model 
specifications include a 5 second advection timestep, Ferrier microphysics, RRTM (Dudhia) log 
wave radiation, Goddard short wave radiation, and Monin-Obukhov surface layer scheme, while 
the North American Regional Reanalysis (NARR; Mesinger 2006) data was used for 
atmospheric initialization and lateral boundary conditions using 3-hourly nudging. Because the 
focus of these simulations are on the sensitivity of PBL-LSM couplings and how they are 
reflected in L-A variables (rather than the absolute accuracy and experiment design), we have 
not included all remaining specifications of LIS-WRF here that are identical for all simulations. 

For each experiment, LIS-Noah and LIS-CLM were run offline (uncoupled) for the 3.5 
year period prior to the start time of IHOP-02 to create equilibrated, or spun-up, land surface 
states for initialization of LIS-WRF. For example, Fig. 2 shows the upper layer (0-10 cm) soil 
water content over the 1-km resolution IHOP-02 domain as generated by the Noah spinup valid 
at 12Z (this procedure was repeated for CLM and for the IHOP-02 June 6 experiments for each 
LSM). The high spatial resolution seen in Fig. 2 is a reflection of the inputs of land cover, 
vegetation, and soil properties available in LIS. Overall, land surface conditions in the ARM- 
SGP region range from highly vegetated and moist in the east to increasingly bare and drier soils 
in the west. 

c. ARM-SGP Observations 

The ARM-SGP program provides surface flux, meteorological, and hydrological 
observations along with atmospheric profiles for a network of sites in and near the winter wheat 
belts of Oklahoma and Kansas. This field experiment has been widely used in previous studies 
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and infonnation on the site locations and characteristics is available at 


www.ann.gov/sites/sgp.stm. 

Radiosondes are launched daily at approximately 1130, 1430, 1730, and 2030 UTC 
(6:30am, 9:30am, 12:30pm, and 3:30pm local solar time) at the SGP central facility at Lamont, 
OK (CF). For this work, radiosonde measurements of temperature, dewpoint, and pressure were 
converted to profiles of G and q at -10 m vertical resolution using standard thennodynamic 
relationships, from which estimates of the height of the PBL (inversion of 6 at the top of the 
mixed layer) were derived. 

The ARM-SGP site employs both Bowen ratio (EBBR) instruments at CF as well as 
numerous extended facilities throughout Kansas and Oklahoma. These data include 30-minute 
average fluxes of net radiation, sensible, latent, and soil heat, along with co-located surface 
radiant temperature, 2-m air temperature, mixing ratio, and wind measurements from 
micrometeorological instrumentation. 

4. Results 

The following sections present mixing diagrams generated from LIS-WRF simulations of 
the IHOP 6 and 12 June 2002 experiments. As described above, generation of these plots and 
derived metrics requires only the diurnal evolution of 0 and q, and mean H s f c and PBLH over the 
time period, which are variables routinely output from LIS-WRF and observed at ARM-SGP. 
a. Mixing diagrams and derived metrics 

The principal controls on the fluxes of heat and moisture from the surface reside in the 
degree of soil moisture and vegetation cover. During spring and summer over the ARM-SGP 
region, there is high spatial variability in each of these (as shown in Fig. 2) that we can use to 
examine LoCo across a range of conditions. Figure 3 presents mixing diagrams from LIS-WRF- 
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Noah and LIS-WRF-CLM simulations on 12 June 2002 for a location with dry soil moisture 
conditions (0.1 1 m 3 m' 3 and 0.08 m 3 m" 3 for Noah and CLM, respectively) in the western and 
bare soil region of the domain. The impacts of coupling the three PBLs to each LSM can be seen 
in the differences in the evolution of 9 and q. 

Qualitatively, the overall shape of the curves indicates little evaporation from the surface 
and significant dry air entrainment into the PBL. The metrics derived from these diagrams are 
also plotted and confirm a high fi s f c and surface heating, while the entrainment flux is primarily 
that of dry air ( f3 ent < -0.50) that results in significant PBL drying and growth throughout the day. 
The large magnitudes of Ah and A/ e confirm that the entrainment fluxes of heat and moisture 
dominate over the surface fluxes, which correspond to rapid and deep PBL growth over a dry 
surface as the values for maximum PBLH are each well over 3 km. All six PBL-LSM 
combinations also indicate a well-established elevated (residual) mixed-layer at the initial time 
that is very unstable, and once this level is reached leads to the explosive and deep PBLs (and 
entrainment fluxes) simulated at this site. In fact, this type of mixing diagram ‘signature’ is 
indicative of an entraimnent feedback loop that supports deep PBL growth, drying of the PBL, 
and desiccation of the surface moisture condition leading to drought if persistent over time (as 
described in Santanello et al. 2007). 

The subtle but significant differences within and between the two diagrams are reflective 
of differences in PBL-LSM coupling. For the Noah LSM (Fig. 3a), the equilibrium created with 
all three PBLs is very similar, with the YSU scheme exhibiting the largest entrainment of heat 
and dry air, and the MYJ scheme the least. This is confirmed in the vertical profile data (not 
shown), where the YSU has an extremely deep maximum PBLH (4.9 km) versus that of the MYJ 
(3.5 km), with the MRF scheme in between (4.2 km). 
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For the CLM simulations (Fig. 3b), there is a noticeable difference in the coupling 
established with the 3 PBLs from that of Noah. All three CLM simulations exhibit greater 
wanning of the PBL (along with comparable drying) than Noah, as reflected in the higher values 
of surface ( /3 S / C > 8.0) and entrainment ((5 ent ~ -0.50) heat fluxes. The sensitivity of CLM to the 
choice of PBL is also greater, as evidenced by the spread between curves in Fig. 3b. The surface 
fluxes produced by CLM are similar for each PBL with very little evaporation taking place, but 
the YSU scheme produces the largest heat and dry air entrainment amounts. It is also important 
to note that the MYJ scheme (as is also the case in Fig. 3a) produces a more slowly growing PBL 
as a result of the initial atmospheric profile being more stable than for the YSU/MRF schemes. 
Overall, the profiles indicate similar maximum PBLH values to the Noah simulations for all 
three PBLs, and are reflective of a drying regime characterized by a desiccated, largely bare soil 
(26 percent vegetation fraction) surface with a deep residual mixed layer that supports significant 
PBL growth and entrainment impacts on the mixed-layer equilibrium. 

Figure 4 presents mixing diagrams for intermediate soil moisture conditions (0.18 m 3 m" 3 
for both Noah and CLM) in the ARM-SGP region as simulated by LIS-WRF using the Noah and 
CLM LSMs. For all PBL-LSM combinations shown in Figs. 4a and 4b, there is a significantly 
different signature of 6 and q evolution than for dry soils. Most significantly, there is little 
diurnal variability in q and a lower dynamic range in 0, which is expected due to the PBL-LSM 
equilibrium created over a more moist and vegetated (40 percent vegetation fraction) surface. 

As described by the metrics, more energy at the surface goes to evaporation (/),/, ) from 
the moister soil at this location, which lowers the amount of surface heating and flux of heat into 
the PBL. There is also less buoyancy and slower PBL growth simulated in each that is reflected 
in the much lower proportion of heat and dry air entrainment (A/,, A/ e ) than for dry soils. The 
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damped evolution of q is a result of the magnitudes of surface evaporation relative to that of dry 
air entrainment, which nearly balance for this location. In this case, there is near zero flux of 
heat into the PBL from entrainment, with some simulations (Noah-YSU and Noah-MRF) 
actually indicating some slightly cooler air mixing through the inversion. Maximum PBLH was 
approximately 1.6 km for the MRF and MYJ simulations, while the YSU was slightly higher 
(~1.8 km; not unexpected given the difference in physics and explicit entraimnent treatment in 
the YSU scheme). Given the close similarity of the three curves and derived fluxes within each 
plot, these results also suggest that the PBL-LSM equilibrium is more significantly impacted by 
the choice of LSM than by the particular PBL scheme employed. However, the impact of each 
LSM and their surface fluxes is diminished at this location due to a strongly stratified initial 
profile that limits PBL growth, and will be discussed further in Section 4c. 

Figure 5 presents mixing diagrams for wet soil conditions (0.32 m 3 m' 3 and 0.40 m 3 m' 3 
for Noah and CLM, respectively) in the eastern portion of the domain, which is also more 
heavily vegetated (> 90 percent vegetation fraction). What is immediately evident for the Noah 
simulations in Fig. 5a is the small range in both 6 and q and the dominance of the moisture flux 
controlled by a nearly freely evaporating surface (low /f/ r ) and limited PBL growth and heat 
entrainment (low /5 ent ). The pattern and fluxes from each of the three PBL schemes is similar, 
and there is very little surface heating and PBL growth (< 1.4 km) for this location as well. One 
interesting feature of this plot is the strong inflection point in q near midday (19 LST); where 
there is a transition from slowly decreasing to rapidly increasing moisture in the PBL. This 
suggests the possibility of a significant horizontal advection component in the system, and will 
be addressed in Section 4e. 
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The CLM simulations (Fig. 5b) produce somewhat different signatures of heat and 
moisture fluxes and PBL evolution at this location. Most notably, surface evaporation is 
dominant, and /3 s / c near zero due to an initial soil condition for CLM that is near saturation. 

There is also a slightly larger diurnal range in temperature than for Noah (higher f5 en t), which 
combined with the near zero surface heat flux produces a higher entrainment (i.e. residual) flux 
estimate. This is partially a result of greater PBL growth in the CLM simulations (PBLH ~ 1.6 
km), though limited as the initial profiles indicate strong stability for both Noah and CLM 
simulations. For all three PBL schemes over the wet surface, A/ e is approximately equal to -1.0 
which indicates the near balance of evaporation with entrainment, but also includes the inflection 
point noted in Fig. 5a. In this case, the diagrams indicate that the choice of LSM is as or more 
important to the simulated L-A equilibrium than the particular PBL scheme employed. This 
follows in that the high moisture availability at the surface and resultant energy balance dominate 
the potential for mixed layer growth and, as a result, the PBL budget over wet surfaces. 
b. PBL budget evaluation 

Mixing diagrams contain a wealth of information that can also be synthesized from an 
energy balance perspective to quantify the variability and accuracy of the processes that govern 
PBL evolution that are difficult to measure. To summarize the L-A processes quantified in these 
mixing diagrams and their sensitivities to different surface moisture conditions, Table 1 lists the 
daily mean heat and moisture fluxes from the surface and entrainment derived from the 
component vectors in Figs. 3-5. The values of these fluxes support the interpretation of the 
diagrams and metrics presented above and highlight the decreasing (increasing) impact of dry air 
entrainment (surface evapotranspiration) from dry to wet soils. Further, these fluxes as a whole 
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define the total heat and moisture budgets of the PBL, and as such define the relative 
contribution of the surface versus atmospheric fluxes. 

Similar to analyses performed for the PILPS experiments (e.g. Chen et al. 1997b, Pitman 
et al. 1999), the relationship of daily mean H and LE can be plotted for each PBL-LSM 
combination. While PILPS was limited to surface fluxes in an offline intercomparison, the 
fluxes derived here also include those of entraimnent, thereby defining the processes contributing 
to the total PBL heat and moisture budget. Following this approach, Fig. 6 shows the 
relationship between sensible and latent heat fluxes from the surface and entraimnent for each of 
the simulations at the dry, intennediate, and wet soil locations shown in Figs. 3-5 and listed in 
Table 1. 

At the dry site, the extremely vigorous PBL development corresponds to the dry air and 
heat entrainment dominating the PBL budget, with is a noticeable sensitivity to the PBL scheme 
employed. The surface fluxes are a much smaller proportion of the total budget, and are similar 
(low LE s f c , high Ilsjc) across all LSMs. For intermediate soils, there is more of a balance between 
surface and entrainment fluxes from which the sensitivities of each to the choice of LSM 
becomes more apparent. However, all six simulations result in a net moisture flux in the PBL of 
approximately -100.0 Wm" 2 , as differences in PBL and surface fluxes cancel out. The sensitivity 
of the PBL budget components to the choice of LSM is more evident for wet soils (Fig. 6c), 
where there is a near balance of evaporation and entrainment ( LE to t ~ 0.0 Wm' 2 ) along with 
limited heating of the PBL (H tot < 200.0 Wm" 2 ), and evidenced by the shallow PBL growth at 
this site. 

Although there are no flux observations at these three locations, the surface fluxes can be 
compared against a theoretical line of constant net radiation (500 Wm' 2 ), which shows how 


22 



available energy at the surface increases for all simulations over wet soils due to changes in 
albedo and surface temperature. This also highlights that the surface available energy varies 
between LSMs, particularly at the intennediate and wet sites, and that the partitioning of fluxes 
for a given LSM varies due to the choice of PBL coupling. Overall, the transition of surface 
fluxes to higher LE s f c and lower H s f c from dry to wet conditions is as expected, and in effect the 
patterns observed in these diagrams describe a transition in the PBL from a drying (entrainment 
dominated) regime to a moistening (surface dominated) regime. 
c. Integrative diagnostics of LoCo 

The vectors and fluxes derived in Figs. 3-6 are, from a slightly broader perspective, 
reflected in two observable properties of the system that are a direct function of the L-A 
equilibrium generated by each PBL-LSM coupling. First, the forcing from the land surface (// v / c 
in Eq. 1) is best represented by the evaporative fraction {EF = LE s f c / ( H s f c + LE s f c )), which is a 
function of the flux of heat and moisture from the land to the atmosphere that contributes to the 
buoyancy and evolution of the PBL. EF is similar to the Bowen ratio but nonnalized for 
incoming available energy, and is sensitive to soil moisture availability as it controls the surface 
flux partitioning. The second integrative property that is sensitive to the PBL-LSM coupling is 
PBL height ( PBLH in Eq. 1), as it is a direct function of the fluxes (most notably heat and dry air 
entrainment) that detennine PBL evolution. 

Combined, the relationship between daily mean EF and maximum PBLH can be thought 
of as describing the amount of surface forcing generated by a LSM versus what the response of 
the coupled PBL scheme is relative to those fluxes. Figure 7 shows an example of this 
relationship as simulated by the LIS-WRF model for the dry, intermediate, and wet soil sites 
depicted in Figs. 3-5. As was shown in Fig. 3, there is significant PBL growth and entrainment 


23 



fluxes over dry soils that are supported by the low EF and very high PBLH values seen here. 

The Noah and CLM simulations with the MYJ PBL showed a slightly different evolution of 0 
and q and lower entrainment rates (A/, and A/ e ) than the YSU/MRF schemes, which is reflected in 
the lower maximum PBLH reached despite having similar low values of EF and high surface 
heating. In this case, the PBL scheme (atmosphere) limits the impact and dampens the forcing 
from the land surface. 

It should be noted that the noticeable differences in the evolution of 0 lq and EFIPBLH 
produced by the YSU/MRF versus the MYJ scheme thus far are not unexpected given the type of 
turbulence closure employed by each scheme. The non-local diffusion parameterized in 
YSU/MRF is known to be more accurate for convective PBL evolution and depth, while the 2.5 
order local closure of the MYJ scheme simulates the nocturnal PBL more accurate and typically 
underestimates the convective PBL (Stensrud 2007). In addition, CLM is traditionally coupled 
with CCSM, which employs a non-local K theory in the PBL similar to that of YSU/MRF. 

For dry and intennediate soil moistures, the CLM simulations produce a slightly lower 
EF (higher H s f c ) than the corresponding Noah runs. This is partly due to the CLM spinup 
initializing a slightly drier soil than Noah for these locations, but also due to the differences in 
LSM physics controlling evaporation. For wet soils, the reverse is true where EF is slightly 
higher in CLM due to slightly higher initial soil water content from the CLM spinup (i.e. despite 
that the spinups used identical atmospheric forcing as for Noah, the model climatology of CLM 
tends to produce higher EF on this date). It is also important to note that the variability in EF 
between Noah and CLM shown here corresponds directly to the differences in V S f C that were 
discussed in Figs. 3-5. 
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For intermediate soils, PBL growth is limited (< 2.0 km) for Noah and CLM despite the 
fact that EF is low and there is strong surface heating (comparable to the dry site) simulated by 
each. This is due to the atmospheric stability and thennal stratification over this site being 
significantly more stable than over the dry site. Specifically, the initial (12Z) stability in the 
lower 3 km of the atmosphere over the intennediate site is approximately 2.6 K/km, while the 
dry site exhibits a much weaker inversion of 4.3 K/km and a deep elevated residual layer. It 
follows (as shown in Fig. 4) that there is only moderate entrainment of heat or moisture, and 
therefore the atmosphere is the dominant control on PBL evolution for the intermediate site. As 
such, PBL evolution is largely insensitive to the choice of LSM, and due to the strong stability is 
also largely insensitive to the choice of PBL scheme and turbulence parameterization employed 
by each. 

For wet soils, there is slightly higher PBLH simulated by CLM despite having a wetter 
surface and higher EF than Noah. This follows with Figs. 5 and 6c, which showed a greater rise 
in CLM mixed layer temperature and entrainment into the PBL than in the Noah simulations. 
There is also slightly more variability in the co-evolution of 0 and q between simulations due to 
the choice of PBL scheme than for intermediate soils. However, the most significant differences 
occur in EF simulated from CLM and Noah, which along with the results of Figs. 5 and 6c 
suggests that the choice of LSM is the more critical component of L-A coupling for this site. As 
such, evaluating the relationship of EF and PBLEl enables the relative strengths and weaknesses 
of the schemes to be identified in terms of observable, integrative properties of the L-A system 
established by each PBL-LSM coupling. There is, as expected, greater sensitivity of EF to the 
choice of LSM and moisture regime, while PBLHv aries more significantly between PBL 
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schemes and particularly for drier soils and when the impacts of the entrainment feedbacks are 
maximized. 

Overall, the series of plots presented in Sections 4a-c have demonstrated the power and 
relative ease of using mixing diagrams to evaluate L-A interactions, as well as their sensitivities 
to differences in PBL-LSM couplings and surface and atmospheric conditions. When combined, 
the three approaches presented above (mixing diagrams, PBL budgets, and EF vs. PBLH) 
provide a comprehensive analysis of the processes governing LoCo in a manner that synthesizes 
the complex interactions and feedbacks into a quantitative and observable framework, while 
offering the ability to directly evaluate the PBL and LSM fluxes and schemes, and sensitivities of 
each. 

d. Mixing diagrams with observations 

This approach can now be supplemented with observations in order to evaluate these 
simulations in the context of the impacts and accuracies of different PBL-LSM couplings. It 
should be noted that it is not the goal of this study to perform an intensive evaluation of the 
physics of the PBL and LSM schemes employed here, but rather to use these experiments to 
demonstrate a framework to further evaluate and understand any coupled modeling system. 

Mixing diagrams from LIS-WRF simulations are presented in Fig. 8 with observations 
made at the ARM-SGP Central Facility (CF; El 3) on 6 June 2002. The Noah simulations with 
all three PBLs generally capture the observed evolution of 0 and q. A closer examination of the 
curves shows that the shift of slope towards negative q (due to dry air entrainment) occurs near 
midday in both the simulations and observations, and there is general agreement in the patterns 
of 6 and q with observations throughout the afternoon including a small increase in q at 23 Z 
when PBL growth and entrainment shuts down. In addition, despite that the metrics show that 
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the Noah model simulates higher (5 s j c values the PBL heights generated by the YSU, MYJ, and 
MRF schemes (1.4, 1.6, and 1.3 km, respectively) are close to observed (1.5 km). This suggests 
that the impact of surface fluxes generated by the LSM is minimal in influencing the evolution of 
PBL properties. 

In contrast, the diurnal patterns of the CLM simulations tend to diverge appreciably from 
the observations and from Noah, and show significant early moistening followed by rapid drying 
due to entrainment. This is a situation where the 12Z/00Z and daily mean values are reasonable, 
but the actual hourly evolution of 0 and q (as well as the mean fluxes) do not reflect reality. 

Once again, all three simulations in Fig. 8b overestimate (5 s f c and underestimate A/„ while the 
YSU scheme produces similar j 6 ent and A/ e values to those observed. The MRF simulation 
produces the deepest PBL (1.6 km), reflected in larger production of dry air entrainment than the 
other two schemes and observations. As was the case for Noah, the slopes of the entrainment 
vectors are quite close to that observed, indicating that regardless of the accuracy of surface flux 
partitioning (by the LSM) or magnitude of entrainment fluxes the ratio of dry air to heat being 
entrained is simulated quite well by all three PBL schemes, as is the maximum PBLH. 

The mixing diagrams for the ARM-SGP facility at Plevna, KS (E4; north of the CF and 
slightly drier) are presented in Fig. 9 for the Noah and CLM simulations. The signatures of 
6 and q are roughly similar to those in Fig. 8, with the Noah simulations noticeably more 
consistent and closer to observations than CLM. From Noah and observations there is less 
moistening of the mixed layer (through evaporation) and consequently higher surface heat fluxes 
than observed (/3 s f c > 2.0), but once again the ratio of dry air being entrained is very close to that 
observed ((3 ent = -0.25) . 
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The CLM simulations show more sensitivity to the PBL scheme, but are closer to 
observed than in Fig. 8b. The YSU and MRF results are initially offset as too warm and moist, 
but the 0 and q evolution throughout the remainder of the day is very close to that observed. The 
MYJ simulation exaggerates the slight morning moistening of the PBL and as a result remains 
too moist throughout the day, which limits PBL growth (< 1.6 km) compared to observations 
(1.8 km). All three PBLs simulate an exaggerated diurnal cycle (particularly in 6 ), and as a 
result also overestimate the residual vector (V en t) and the proportion of heat versus dry air being 
entrained (Pent » -0.25) relative to that observed. 

The PBL budget and EF/PBLH diagnostics can now be examined to further understand 
the nature of the L-A coupling created by each of the PBL-LSM pairings and how they compare 
to observations at these two sites. Fluxes comprising the PBL budget of heat and moisture from 
the E13 and E4 sites (shown in Figs. 8 and 9) are shown in Fig. 10. Also included are the 
observed fluxes derived from the mixing diagrams and lines of constant energy (H s f c + LE s f c ) 
observed at the land surface. At the E13 site, the surface fluxes from each simulation align 
below the observed available energy, with the flux partitioning biased (~ 100 Win" 2 ) towards 
higher H s f c and lower LE s f c than observed (as reflected in the slopes of V S f C in Fig. 8). On the 
other hand, the entrainment and total PBL fluxes indicate a greater spread in energy partitioning, 
with an underestimation of H ent that is primarily greater than the overestimation of heat flux from 
the surface. When compared against the observed total fluxes (i.e. budget), this results in four of 
the six simulations underestimating the net total heat flux into the PBL at this site. Likewise, 
there is more dry air entrainment than evaporation of moisture into the PBL, leading to a 
negative input of moisture into the PBL from all six PBL-LSM combinations. 
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At the E4 site (Fig. 10b) a similar pattern emerges in that the available energy and 
evaporation at the surface are underestimated in all simulations, but unlike E13 there is a clear 
sensitivity to the choice of LSM. This follows for entrainment as well, where Noah (and in 
particular Noah-YSU) simulates dry and heat entraimnent fluxes very close to those observed. 

As the mixing diagrams in Fig. 9 indicate, CLM produces too much heat entrainment at the 
expense of dry air. Therefore, because all six PBL-LSM couplings produce too little 
evaporation, the total moisture budget from Noah (CLM) is too dry (wet), while the heat input to 
the PBL is slightly overestimated in each. 

With the addition of observations to the EF/PBLH diagnostics, it can also now be 
ascertained if each PBL-LSM coupling produces the ‘correct’ answer (e.g. PBLH) despite flaws 
in the representation of specific L-A processes (e.g. surface fluxes, entrainment). Figure 1 1 
presents the relationship of EF and maximum PBLFl at the E13 and E4 sites along with 
observations. As suggested by the PBL budget results at E13, there is not a clear indication that 
the PBL or LSM choice is more important than the other, and the PBL growth is sensitive to the 
precise nature of the PBL-LSM coupling. When combined with Fig. 11a, however, it can be 
ascertained that the two outlier points in Fig. 10a that overestimate the total heat budget (~ 400 
WnT 2 ) are also the same PBL-LSM couplings (Noah-MYJ and CLM-MRF) that are outliers in 
terms of EF and PBLF[. The remaining four simulations produce less PBL growth, and all 
simulations estimate significantly lower surface evaporation than observed. 

The relationships at the E4 site show a clear sensitivity to the choice of LSM, with CLM 
producing higher evaporation (in response to a slightly wetter soil than Noah) throughout the day 
and in return diminished PBL growth compared to Noah and observations. This is confirmed in 
the PBL budget and mixing diagram results, where CLM underestimates the drying due to 
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entrainment. Once again all PBL-LSM couplings underestimate the EF, and at this site the 
observed maximum PBLH is not reached by any of the simulations. This may be due to the fact 
that the surface available energy is too low, thereby limiting evaporation and also the buoyant 
energy at the surface. 

Overall, the combination of the three approaches (Figs. 8-11), with their foundation in the 
mixing diagram approach and the addition of observations, can be useful in diagnosing the 
accuracy and sensitivity of land surface and PBL state variables, fluxes, and bulk PBL 
properties. For example, while CLM does a poorer job than Noah of simulating the diurnal 
evolution of 0 and q at E13 and E4, there is clearly more sensitivity to the choice of LSM at site 
E4 (as evidenced in Figs. 10 and 1 1). Furthermore, flux biases (such as that shown in Fig. 10a 
for H ent ) can be diagnosed that indicate errors in the amount of energy supplied at the land 
surface and/or into the PBL through entrainment created by a particular PBL-LSM coupling. 
These are all important steps towards greater and complete understanding and quantification of 
the components of LoCo. 

5. Discussion 

The mixing diagram theory presented by Betts (1992) also supports inclusion of a 
horizontal advection vector. As many studies have shown (Kustas and Brutsaert 1987; Peters- 
Lidard and Davis 2000; Santanello et al. 2005), one of the main limiting factors (other than 
entrainment) in closing the heat and moisture budgets of the PBL is advection. Here, horizontal 
advection of heat and moisture can be represented by a vector (V a dv) in the same manner as the 
surface and entrainment fluxes, and represents the horizontal flux of heat and moisture over the 
time period. The contribution of advection is calculated and then added to the surface flux 
vector (Vadv + V S fc), with the new residual representing the entrainment flux as depicted in Fig. 
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12. As such, the surface flux vector is unaffected by the addition of advection, while 
entrainment clearly is impacted to a degree detennined by the magnitude and direction of the 
advection vector. The advection Bowen ratio (/L/,,) and component fluxes (H ac / V , LE a dv ) are then 
computed in analogous fashion to the surface and entrainment contributions to the PBL budget. 
The high-resolution of LIS-WRF output (1 km; hourly) makes it relatively easy to calculate 
hourly advection estimates using a finite differencing approach and to re-generate mixing 
diagrams that include all three PBL budget components. 

Figure 13 presents the mixing diagrams and for the wet soil site (as shown in Fig. 5) after 
the inclusion of the advection vectors. In the Noah simulations, the advection is in the cold and 
moist direction but rather small relative to the surface and entrainment vectors. In contrast, the 
CLM simulations produce more significant advection fluxes that are opposite in sign (wann and 
dry) to those from Noah. As such, the CLM advection is acting in the same direction as 
entrainment thereby lowering the magnitude of the residual vector (Vent). This becomes evident 
in the flux values and PBL budgets plotted in Fig. 14. The original CLM mixing diagrams 
(ignoring advection) and results from Figs. 5-7 indicate rather high estimates of entrainment 
fluxes considering the high surface evaporation rate which limited surface heating and PBL 
growth at this site. In the new diagram, it is evident that advection is contributing to the 
warming and drying of the mixed layer, and as a result the surface evaporation can be thought of 
as being nearly balanced by entrainment and advection (i.e. LE tot ~ 0.0; Ai eient + A/ eia( / v ~ -1.0). 

The magnitude of the advection vector for the CLM makes it an important addition to the 
mixing diagram approach and interpretation of its derived metrics and fluxes. However, this is 
not the case for Noah or the other sites and results presented earlier, where advection and the 
resultant impacts on L-A coupling processes and PBL budgets were minimal. While the 
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inclusion of advection does not impact the surface or total energy budget value, or the 
relationship between EF and PBLH, it should still be considered on a case by case basis as it may 
comprise a large component of the residual vector and can easily be applied as shown here. 

In terms of the residual (entrainment) vector itself, it can best be put in terms of an 
‘atmospheric response’ vector in that it represents the full sum of atmospheric fluxes and 
contribution to PBL evolution. While entrainment was shown to be dominant for the case 
studies selected here, this list includes processes such as advection, radiative flux divergence, 
compressional warming, and moist processes such as condensation/evaporation at the top of the 
PBL. The clear sky focus and approach taken here can therefore be expanded (as for advection) 
to account for these processes as additional components of the residual vector, should they be 
significant. 

An important application of mixing diagrams is the ability to evaluate coupled systems in 
terms of the accuracy of their component surface and PBL schemes. In comparing the evolution 
of 0 and q, fluxes, and bulk PBL properties (e.g. PBLH) to other models and observations, 
deficiencies in the various schemes’ production of turbulent fluxes can be pinpointed. Typically, 
LSMs are evaluated offline and against individual variables or fluxes. Using the approach 
presented here, however, modifications to LSM parameters or physics can be implemented and 
their impacts on simulated surface variables and fluxes can be evaluated along with the 
corresponding atmospheric response in the fonn of PBL evolution and budgets (i.e. the residual 
vector). Using the E4 site above as an example, modifications to increase surface available 
energy and evaporation in the LSM scheme would certainly impact the turbulent fluxes and 
evolution of the PBL (including non-linear feedback effects), resulting in a new equilibrium 
between L-A processes as reflected in the diagrams. The impact and sensitivity of each coupled 
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system to modifications to the LSM (or PBL) can therefore be quantified in tenns of changes to 
the vector components and PBL budgets. 

The mixing diagram approach can also be used to identify poorly understood L-A 
feedback regimes in coupled systems. For example, the feedback of dry air entrainment on the 
evolution of surface fluxes can be a significant detenninant of the diurnal L-A equilibrium 
created between the PBL and land surface (Santanello et al. 2007). Figures 3-7 showed the 
differences in this equilibrium across surface moisture regimes, but also highlight the feedback 
of dry air entrainment in a) maintaining a nearly constant evaporative demand and moisture 
budget in the PBL at the intermediate and wet sites, and b) supporting extremely vigorous PBL 
growth and dry air entrainment due to the presence of an elevated residual mixed layer at the dry 
site. These feedbacks tend to be self-perpetuating, so identification of each is important in 
evaluating various PBL-LSM couplings and their limitations. Further, the stratification over 
surface moisture regimes could easily be perfonned over varying land cover, vegetation, soil, or 
atmospheric properties to evaluate the sensitivities of the coupled system to a wide range of 
conditions. 

There are potential issues relating to scale when applying these methods that must be 
addressed as well. The LIS-WRF experiments were run at very high spatial resolution (1 km), 
and for the analyses presented here a single grid cell, or column, was pulled out nearest to the 
site of interest. Clearly, there is heterogeneity in land and atmospheric properties around each 
cell that may impact the evolution of PBL properties and fluxes. Therefore, additional analyses 
were performed using mixing diagrams created from columns taken from up to 25 km 
surrounding each point of interest. It was found that, while there is variability in 0 , q, and PBL 
evolution around each site, differences were small relative to the overall patterns observed and 
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characteristics of the central location, and therefore did not impact the demonstration of the 
approach or interpretation of results. In addition, the surface flux, temperature, and humidity 
observations are point measurements, and best evaluated at that point rather than an areal 
average. Studies and locations that are more heterogeneous may require such an averaging 
approach, however, particularly with regards to PBL properties (e.g. PBLH) that integrate over a 
much larger region. 

6. Conclusions 

The framework and results presented here provide a comprehensive methodology to 
quantify and evaluate the critical processes controlling local L-A coupling. The ability to 
evaluate the full PBL heat and moisture budgets and their flux components is critical to 
identifying the dominant processes involved in LoCo as well as deficiencies in PBL and LSM 
schemes and their interactions and feedbacks. As was shown here, all the information necessary 
for such an analysis is contained in mixing diagrams, wherein the diurnal co-evolution of G and q 
is integrative of the processes controlling PBL growth and the resultant L-A equilibrium 
established. Overall, the combination of mixing diagrams with their derived metrics (e.g. 
entrainment and Bowen ratios), PBL budgets, and integrative diagnostics ( EF vs. PBLH) with 
observations supplies a consistent and practical framework from which to evaluate coupled 
models and parameterization schemes on diurnal time scales. 

The IHOP-02 experiments and sites focused on here cover a wide range of surface 
moisture and atmospheric conditions from which to test the mixing diagram approach and its 
sensitivity to significant variability in PBL and land surface fluxes. The different ‘signature’ of 
each mixing diagram is quite evident, and reflects a strong sensitivity of the derived fluxes and 
metrics to these conditions and the various combinations of PBL-LSMs, and can be evaluated in 
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combination with the PBL budgets and EF/PBLH relationships. For dry soils, the PBL evolution 
and structure is more significantly impacted by the choice of PBL scheme than by the particular 
LSM employed. For intermediate soils, the diagrams indicate that neither the choice of LSM or 
PBL dominates the simulated L-A equilibrium due to a strongly stable atmosphere that 
constrains the system. For wet soils, there is more influence of the choice of LSM on the 
resultant PBL evolution and fluxes. 

The ultimate utility of this approach is in evaluating coupled models and their scheme 
components against observations. The two ARM-SGP sites presented here (El 3 and E4) show 
significant variability in 9/q signatures and derived fluxes due to the specific PBL-LSM coupling 
employed. The evaluation of PBL budget components as well as PBL height suggests that there 
is not enough available energy to support evaporation at the surface and that PBL turbulence is 
underestimated, the degree to which depends on the scheme choice. Modifications to these 
schemes therefore will result in a new L-A equilibrium that can similarly be evaluated against 
the observations. This is a major advantage of the methodology, in providing a truly coupled 
evaluation of the processes involved in terms of observable, integrative properties of the system, 
rather than a traditional one-at-a-time or offline approach. 

The IHOP-02 experiments were provided as an example of how to apply the mixing 
diagram approach to model output and observations. While the focus is on diurnal and local 
scales for convective PBLs, this framework can be easily applied to any coupled model, scales, 
and conditions of interest. Although the focus was on cloud-free and smooth diurnal cycles, the 
mixing diagram framework used in this study includes the full set of governing L-A interactions 
and processes that allows for other applications. For example, the entrainment Bowen ratio 
(p en t), shown to be an important determinant of convective initiation (Betts and Ball 1994) yet 
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difficult to measure, is one of the principal metrics derived from this approach. Such metrics 
would be valuable to understanding the generation of convection in coupled models by 
quantifying the L-A processes and feedbacks that are typically difficult to interpret (Trier et al. 
2004; Holt et al. 2006). 

Ongoing work on advancing the cause of LoCo includes a number of detailed 
experiments and analyses based on the mixing diagram approach. For example, an evaluation of 
different methods to spinup Noah and CLM and initialize LIS-WRF is being performed with 
varying degrees of input forcing and parameter data quality to yield insight into the sensitivity 
and accuracy of various PBL-LSM couplings to the initial conditions. Another ongoing 
experiment is an extended (~ 7 day) regional simulation that will enable LIS-WRF to evolve 
over time from synoptically forced, clear-sky conditions with a drydown period ending with the 
convectively-active 12 June case. As such, a longer-term transition can be seen in the mixing 
diagrams from each day that reflects the changing surface and atmospheric conditions through 
the evolution of 0 and q. 

In addition, this framework will be included in an upcoming GEWEX-GLASS directed 
community-wide set of pilot experiments, where LIS-WRF will serve as the testbed. This study 
will evaluate a large set of coupling diagnostics in order to develop a hierarchical list of coupling 
coefficients for LoCo. These diagnostics include some of the efforts mentioned in Section 2, and 
will cover the range of local scale interactions from surface-PBL (e.g. mixing diagrams) to moist 
processes and convective triggering. LIS-WRF has recently added additional LSM options and 
WRF (v.3) has added another PBL scheme, both of which will be employed in this study along 
with a stand-alone single-column model testbed. 
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Quantification of L- A interactions is particularly important for land surface data 
assimilation and model calibration efforts. While these are relatively young topic of research, 
high-quality remote sensing data (e.g. surface temperature, snow, and soil moisture) can be 
assimilated into LSMs using a variety of techniques and used to calibrate both land and/or 
atmospheric parameters. However, the impacts of these techniques are vastly different for 
offline and coupled models (e.g. Liu et al. 2003, 2004, 2005) due to the addition of L-A 
interactions and feedbacks in the latter. The mixing diagram approach can therefore be an 
important tool in detennining the potential improvement and model sensitivity to assimilation 
and calibration strategies going forward. 

Finally, the greater applicability of this methodology to the LoCo community is not 
limited to modeling studies alone. Recent advances in satellite remote sensing will continue to 
improve the retrieval of PBL and land surface data for a number of applications with global 
coverage and high temporal resolution. This includes the diurnal evolution (due to multiple 
sensors) of variables such as temperature and humidity (MODIS, AIRS), soil moisture (AMSR, 
SMOS), evaporation (MODIS, AIRS), and PBL height (AIRS, CALIPSO). As a result, the 
ability of satellite remote sensing to monitor the PBL and estimate L-A properties and conditions 
will continue to be improved and can be incorporated to the mixing diagram approach to provide 
insight into LoCo across the globe. 

Acknowledgements. This work was supported by the NASA Energy and Water Cycle Study 
(NEWS; PI: Peters-Lidard) and we would like to thank NEWS, ESSIC, and GSFC for helping to 
make the completion of this work possible. In particular, Jim Geiger and Joe Eastman were 
instrumental in providing feedback and activities related to LIS-WRF. We also appreciate the 


37 



past and ongoing collaboration with the LoCo community that has stimulated this work, in 
particular Bert Holtslag, Bart van den Hurk, Paul Houser, and Dara Entekhabi. 


38 



6. References 


Angevine W. M., 1999: F.n train merit results including advection and case studies from the 
Flatland boundary layer experiments. J. Geophys. Res., 104, 30947-30963. 

Anthes, R.A., and T.T. Warner, 1978: Development of Hydrodynamic Models Suitable for Air 
Pollution and Other Mesometerological Studies. Mon. Wea. Rev., 106, 1045-1078. 

Barros A. P., and W. Hwu, 2002: A study of land-atmosphere interactions during summertime 
rainfall using a mesoscale model. J. Geophys. Res., 107, 4227, 
doi:10.1029/2000JD000254. 

Berbery, E., Y. Luo, K. Mitchell, and A. Betts, 2003: Eta model estimated land surface processes 
and the hydrological cycle of the Mississippi Basin. J. Geophys. Res., 108, D22, 8852, 
doi: 10. 1029/2002JD003 192. 

Best, M.J., A. Beljaars, J. Polcher, and P. Viterbo, 2004: A Proposed Structure for Coupling 
Tiled Surfaces with the Planetary Boundary Layer. J. Hydrometeor., 5, 1271-1278. 

Betts, A.K., 1984: Boundary Layer Thermodynamics of a High Plains Severe Storm. Mon. Wea. 
Rev., 112,2199-2211. 

, 1992: FIFE atmospheric boundary layer budget methods. J. Geophys. Res., 97, 1 8523— 

18532. 

, 2000: Idealized model for equilibrium boundary layer over land. J. Hydrometeor., 1, 507- 

523. 

, and J. H. Ball, 1994: Budget analysis of FIFE-1987 sonde data. J. Geophys. Res., 99, 

3655-3666. 

, , A. C. M. Beljaars, M. J. Miller, and P.A. Viterbo, 1996: The land surface- 

atmosphere interaction: A review based on observational and global modeling 


39 



perspectives. J. Geophys. Res., 101, 7209-7226. 

, and P. Viterbo, 2005: Land-surface, boundary layer, and cloud- field coupling over the 

southwestern Amazon in ERA-40. J. Geophys. Res., 110, D14108, 
doi: 10. 1029/2004 JD005702. 

Bosilovich, M. G., J. D. Radakovich, A. da Silva, R. Todling, and F. Verter, 2007: Skin 

temperature analysis and bias correction in a coupled land-atmosphere data assimilation 
system. J. Met. Soc. Japan, 85A, 205-228. 

Brubaker, K. L., and D. Entekhabi, 1996: Analysis of feedback mechanisms in land-atmosphere 
interaction. Water Resour. Res., 32, 1343-1357. 

Chen, F., Z. Janie, and K. Mitchell, 1997a: Impact of atmospheric-surface layer 

parameterizations in the new land-surface scheme of the NCEP mesoscale Eta numerical 
model. Bound. -Layer Meteor., 85, 391-421. 

Chen, T.H., A. Henderson-Sellers, P.C.D. Milly, A.J. Pitman, A.C.M. Beljaars, J. Polcher, F. 

Abramopoulos, A. Boone, S. Chang, F. Chen, Y. Dai, C.E. Desborough, R.E. Dickinson, 
L. Dumenil, M. Ek, J.R. Garratt, N. Gedney, Y.M. Gusev, J. Kim, R. Koster, E.A. 
Kowalczyk, K. Laval, J. Lean, D. Lettenmaier, X. Liang, J.F. Mahfouf, H.T. 
Mengelkamp, K. Mitchell, O.N. Nasonova, J. Noilhan, A. Robock, C. Rosenzweig, J. 
Schaake, C.A. Schlosser, J.P. Schulz, Y. Shao, A.B. Shmakin, D.L. Verseghy, P. Wetzel, 
E.F. Wood, Y. Xue, Z.L. Yang, and Q. Zeng, 1997b: Cabauw Experimental Results from 
the Project for Intercomparison of Land-Surface Parameterization Schemes. J. Climate, 
10, 1194-1215. 

Cheng, W.Y.Y., and W.J. Steenburgh, 2005: Evaluation of Surface Sensible Weather Forecasts 
by the WRF and the Eta Models over the Western United States. Wea. Forecasting, 20, 


40 



812-821. 


Cleugh H. A., M. R. Raupach, P. R. Briggs, and P. A. Coppin, 2004: Regional-scale heat and 

water vapour fluxes in an agricultural landscape: An evaluation of CBL budget methods 
at OASIS. Bound. -Layer Meteorol., 110,99-137. 

Dai Y. J., Coauthors,, 2003: The Common Land Model. Bull. Amer. Meteor. Soc, 84, 1013— 
1023. 

Desai, A. R., Davis, K. J., Senff, C. J., Ismail, S., Browell, E. V., Stauffer, D. R., and 

Reen, B. P.: 2006: A Case Study on the Effects of Heterogeneous Soil Moisture on 
Mesoscale Boundary Layer Structure in the Southern Great Plains, U.S.A. Part I: Simple 
Prognostic Model. Bound. -Layer Meteorol., 119, DOI: 10.1007/sl0546-005-9024-6. 

Diak G. R., 1990: Evaluation of heat flux, moisture flux and aerodynamic roughness at the land 
surface from knowledge of the PBL height and satellite derived skin temperatures. Agric. 
For. Meteor., 22, 505-508. 

Dirmeyer, P. A., A. J. Dolman, and N. Sato, 1999: The Global Soil Wetness Project: A pilot 
project for global land surface modeling and validation. Bull. Amer. Meteor. Soc., 80, 
851-878. 

Dirmeyer P. A., Z. Guo, and X. Gao, 2004: Comparison, validation, and transferability of eight 
multiyear global soil wetness products. J. Hydrometeor., 5, 1011-1033. 

Dolman A., J. Gash, J. Goutorbe, Y. Kerr, T. Lebel, S. Prince, and J. Strieker, 1997: The role of 
the land surface in Sahelian climate: HAPEX-Sahel results and future research needs. J. 
Hydro!, 188/189, 1067-1079. 

Ek M. B., and A. A. M. Holtslag, 2004: Influence of soil moisture on boundary layer cloud 
development. J. Hydrometeor., 5, 86-99. 


41 



Ek M. B., K. E. Mitchell, Y. Lin, E. Rogers, P. Grunmann, V. Koren, G. Gayno, and J. D. 

Tarpley, 2003: Implementation of Noah land surface model advances in the National 
Centers for Environmental Prediction operational mesoscale Eta Model. J. Geophys. Res., 
108, 8851, doi: 10.1029/2002JD003296. 

Eltahir E. A., 1998: A soil moisture-rainfall feedback mechanism: 1. Theory and observations. 
Water Resour. Res., 34, 765-776. 

Entekhabi, D., and K. L. Brubaker, 1995: An analytic approach to modeling land-atmosphere 
interaction. 2: Stochastic extension. Water Resour. Res., 31, 633-643. 

Entekhabi D., Coauthors, 1999: An agenda for land surface hydrology research and a call for the 
second international hydrology decade. Bull. Amer. Meteor. Soc., 80, 2043-2058. 

Findell K. L., and E. A. Eltahir, 2003a: Atmospheric controls on soil moisture-boundary layer 
interactions: Three-dimensional wind effects. J. Geophys. Res., 108, 8385, 
doi: 10. 1029/2001 JD0015 15. 

Findell K. L., and E. A. Eltahir, 2003b: Atmospheric controls on soil moisture-boundary layer 
interactions. Part I: Framework development. J. Hydrometeor., 4, 552-569. 

Garratt, J. R., 1992: The Atmospheric Boundary Layer. Cambridge University Press, 316 pp.. 

Gu, L., T. Meyers, S. G. Pallardy, P. J. Hanson, B. Yang, M. Heuer, K. P. Hosman, J. S. Riggs, 

D. Sluss, and S. D. Wullschleger, 2006: Direct and indirect effects of atmospheric 
conditions and soil moisture on surface energy partitioning revealed by a prolonged 
drought at a temperate forest site, J. Geophys. Res., 1 1 1, D 16 102, 
doi: 1 0. 1 029/2006 JD007 161. 

Hacker, J.P., and D. Rostkier-Edelstein, 2007: PBL State Estimation with Surface Observations, 
a Column Model, and an Ensemble Filter. Mon. Wea. Rev., 135, 2958-2972. 


42 



Hacker, J.P., and C. Snyder, 2005: Ensemble Kalman Filter Assimilation of Fixed Screen-Height 
Observations in a Parameterized PBF. Mon. Wea. Rev., 133 , 3260-3275. 

Henderson-Sellers, A., Z.F. Yang, and R. Dickinson, 1993: The Project for Intercomparison of 
Fand-surface Parameterization Schemes. Bull. Amer. Meteor. Soc., 74 , 1335-1349. 

Hess R., 2001: Assimilation of screen-level observations by variational soil moisture. Meteor. 
Atmos. Phys., 77, 145-154. 

Hogue, T. S., F. Bastidas, H. Gupta, S. Sorooshian, K. Mitchell, and W. Emmerich, 2005: 
Evaluation and Transferability of the Noah Fand Surface Model in Semiarid 
Environments. J. Hydrometeorol., 6 , 68-84. 

Holt, T. R., D. Niyogi, F. Chen, K. Manning, and M. A. FeMone, and A. Qureshi, 2006: Effect 
of land-atmosphere interactions on the IHOP 24-25 May 2002 convection case. Mon. 
Wea. Rev., 134 , 113-133. 

Holtslag, A. A. M., G. J. Steeneveld, and B. J. H. van de Wiel, 2007: Role of land-surface 
temperature feedback on model performance for the stable boundary layer. 

Bound. -Layer Meteorol, 125 , 361-376. 

Hong, S.Y., and H.F. Pan, 1996: Nonlocal Boundary Fayer Vertical Diffusion in a Medium- 
Range Forecast Model. Mon. Wea. Rev., 124 , 2322-2339. 

Hong, S.Y., Y. Noh, and J. Dudhia, 2006: A New Vertical Diffusion Package with an Explicit 
Treatment of Entrainment Processes. Mon. Wea. Rev., 134 , 2318-2341. 

Jacobs, C., and H. De Bruin, 1992: The Sensitivity of Regional Transpiration to Fand-Surface 
Characteristics: Significance of Feedback. J. Climate, 5 , 683-698. 

Janjic, Z. I., 2001: Nonsingular implementation of the Mellor-Yamada level 2.5 scheme in the 

NCEP meso model. Technical Report 437, National Centers for Environmental Prediction 
Office. 


43 



Kim C. P., and D. Entekhabi, 1998: Feedbacks in the land-surface and mixed-layer energy 
budgets. Bound. -Layer Meteor., 88, 1-21. 

Koster R. D., P. A. Dinneyer, A. N. Hahmann, R. Ijpelaar, L. Tyahla, P. Cox, and M. J. Suarez, 
2002: Comparing the degree of land-atmosphere interaction in four atmospheric general 
circulation models. J. Hydrometeor, 3 , 363-375. 

Koster R. D., Coauthors, 2004: Regions of strong coupling between soil moisture and 
precipitation. Nature, 306 , 1 138-1 140. 

Kumar, S. V., C. D. Peters-Lidard, J. L. Eastman, and W.-K. Tao, 2008: An integrated high 
resolution hydrometeorological modeling testbed using LIS and WRF. Environmental 
Modelling and Software, 23 , 169-181. 

Kustas W. P., and W. Brutsaert, 1987: Virtual heat entrainment in the mixed layer over very 
rough terrain. Bound. -Layer Meteor., 38, 141-157. 

Lawrence, D.M., and J.M. Slingo, 2005: Weak Land- At mo sphere Coupling Strength in 
HadAM3: The Role of Soil Moisture Variability. J. Hydrometeor., 6, 670-680. 

Liu, Y., L. A. Bastidas, H. V. Gupta, and S. Sorooshian, 2003: Impacts of a Parameterization 

Deficiency on Offline and Coupled Land Surface Model Simulations. J. Hydrometeorol., 
4 , 901-914. 

Liu, Y., H. V. Gupta, S. Sorooshian, L. A. Bastidas, and W. J. Shuttleworth, 2004: Exploring 
Parameter Sensitivities of the Land Surface Using a Locally Coupled Land-Atmosphere 
Model. J. Geophys. Res., 109 , 21101-21114. 

Liu, Y., H. V. Gupta, S. Sorooshian, L. A. Bastidas, and W. J. Shuttleworth, 2005: Constraining 
Land Surface and Atmospheric Parameters of a Locally Coupled Model Using 
Observational Data. J. Hydrometeorol., 6, 156-172. 


44 



Margulis, S. A., and D. Entekhabi, 2001: Feedback between the Land Surface Energy Balance 
and Atmospheric Boundary Layer Diagnosed through a Model and its Adjoint. J. 
Hydrometeorol., 2 , 599-620. 

Medeiros, B., A. Hall, and B. Stevens, 2005: What Controls the Mean Depth of the PBL? J. 
Climate, 18 , 3157-3172. 

Mellor, G. L., and T. Yamada, 1982: Development of a turbulence closure model for geophysical 
fluid problems. Rev. Geophys. Space Phys., 20 , 851-875. 

Mesinger, F., G. DiMego, E. Kalnay, K. Mitchell, P.C. Shafran, W. Ebisuzaki, D. Jovic, J. 

Woollen, E. Rogers, E.H. Berbery, M.B. Ek, Y. Fan, R. Grumbine, W. Higgins, H. Li, Y. 
Lin, G. Manikin, D. Parrish, and W. Shi, 2006: North American Regional Reanalysis. 
Bull. Amer. Meteor. Soc., 87, 343-360. 

Michalakes, J., S. Chen, J. Dudhia, L. Hart, J. Klemp, J. Middlecoff, and W. Skamarock, 2001: 
Development of a next generation regional weather research and forecast model. 
Developments in Telecomputing: Proceedings of the Ninth ECMWF Workshop on the 
use of high performance computing in meteorology. Singapore, pp. 269-276. 

Molod, A., H. Salrnun, and D.W. Waugh, 2004: The Impact on a GCM Climate of an Extended 
Mosaic Technique for the Land- At mo sphere Coupling. J. Climate, 17 , 3877-389 1. 

Oke T. R., 1987: Boundary Layer Climates. 2d ed. Routledge, 435 pp. 

Pan, H.-L., and L. Mahrt, 1987: Interaction between soil hydrology and boundary-layer 
development. Bound. -Layer Meteor., 38, 185-202.. 

Peters-Lidard C., and L. H. Davis, 2000: Regional flux estimation in a convective boundary layer 
using a conservation approach. J. Hydrometeor, 1 , 170-182. 

Pitman A. J., Coauthors, 1999: Key results and implications from phase 1(c) of the Project for 


45 



Intercomparison of Land-surface Parameterization Schemes. Climate Dyn., 15, 673-684. 

Polcher J., Coauthors, 1998: A proposal for a general interface between land-surface schemes 
and general circulation models. Global Planet. Change, 19, 261-276. 

Reen, B. P., D. R. Stauffer, K. J. David, and A. Desai, 2006: A Case Study on the Effects of 
Heterogeneous Soil Moisture on Mesoscale Boundary Layer Structure in the Southern 
Great Plains, U.S.A. Part II: Mesoscale Modeling. Bound. -Layer Meteorol., 120, DOI: 
10. 1007/s 10546-006-9056-6. 

Reichle, R. H., R. D. Koster, P. Liu, S. P. P. Mahanama, E. G. Njoku, and M. Owe, 2007: 
Comparison and assimilation of global soil moisture retrievals from the Advanced 
Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) and the 
Scanning Multichannel Microwave Radiometer (SMMR). J. Geophys. Res., 112, 
D09108. 

Rhodin A., F. Kucharski, U. Callies, D. P. Eppel, and W. Wergen, 1999: Variational analysis of 
effective soil moisture from screen-level atmospheric parameters: Application to a short- 
range weather forecast model. Quart. J. Roy. Meteor. Soc., 125, 2427-2448. 

Rodell, M. and P. R. Houser, 2004: Updating a land surface model with MODIS -derived snow 
cover. Journal of Hydrometeorology , 5, 1064-1075. 

Santanello, J. A., M. A. Friedl, and M. Ek, 2007: Convective Planetary Boundary Layer 

Interactions with the Land Surface at Diurnal Time Scales: Diagnostics and Feedbacks. 
J. Hydrometeor., 8, 1082-1097. 

Santanello, J. A., M. A. Friedl, and W. P. Kustas, 2005: Empirical Investigation of Convective 
Planetary Boundary Layer Evolution and its Relationship with Land Surface Properties 
and Processes. J Appl. Meteorol., 44, 917-932. 


46 



Seuffert G., H. Wilker, P. Viterbo, M. Drusch, and J. F. Mahfouf, 2004: The usage of screen- 
level parameters and microwave brightness temperature for soil moisture analysis. J. 
Hydrometeor, 5 , 516-531.Sorbjan 1995 

Sorbjan, Z., 1995: Toward evaluation of heat fluxes in the convective boundary layer. J. Appl. 
Meteor ol., 34 , 1092-1098. 

Steeneveld, G.J., B.J.H. van de Wiel, and A.A.M. Holtslag, 2006: Modeling the Evolution of the 
Atmospheric Boundary Layer Coupled to the Land Surface for Three Contrasting Nights 
in CASES-99. J. Atmos. Set., 63 , 920-935. 

Stensrud, D. J., 2007: Parameterization Schemes: Keys to Understanding Numerical Weather 
Prediction Models. Cambridge University Press, 459 pp. 

Stull, R. B., 1988: An Introduction to Boundary Layer Meteorology. Kluwer Academic, 666 pp. 

Trier, S. B., F. Chen, and K. W. Manning, 2004: A study of convection initiation in a mesoscale 
model using high-resolution land surface initial conditions. Mon. Wea. Rev., 132 , 2954- 
2976. 

Troen, I., and L. Mahrt, 1986: A simple model of the atmospheric boundary layer: Sensitivity to 
surface evaporation. Bound. -Layer Meteor., 37, 129-148.. 

van den Hurk, B., A. Holtslag, and C. D. Peters-Lidard, 20005: GLASS and GABLS workshop 
on local land-atmosphere coupling - 19/21 Sep. 2005: Workshop Summary and LoCo 
Implementation Plan. 

Weckworth T. M., Coauthors, 2004: An overview of the International FLO Project (IHOP 2002) 
and some preliminary highlights. Bull. Amer. Meteor. Soc, 85 , 253-277. 


47 



List of Figures 


Figure 1: a) Diurnal evolution (12-00Z; solid line) of 2m-potential temperature (C p 6) vs. 2m- 
specific humidity (Lq) from a representative day during June 2002 in the U.S. Southern Great 
Plains as simulated by a coupled mesoscale model. The annotations on the plots depict the 
vector component contributions of surface and entrainment fluxes that are obtained using the 
mixing diagram approach, b) 12-hour time series of 2m-potential temperature (0) and 2m- 
specific humidity (q) corresponding to Fig. la. c) 12-hour time series of PBL height ( PBLH ) 
and surface sensible heat flux (H s j c ) which are used in calculating the vector and flux components 
in Fig. la. 

Figure 2: Soil moisture (volumetric: m 3 m" 3 ) in the upper (0-10cm) layer valid at 12Z on 12 
June 2002 as simulated from a 2.5 year spinup of the Noah model over the 1-km LIS-WRF 
domain in the SGP for the IHOP-02 experiment. The ARM-SGP Central Facility (CF) at 
Lamont, OK is also shown along with the dry, intermediate, and wet soil locations presented in 
Figs. 3-5. 

Figure 3: Diurnal co-evolution (12-00Z) of 2m-specific humidity (Lq) and 2m-potential 

temperature ( C p 6 ) on 12 June 2002 as simulated by LIS-WRF during dry soil moisture 
conditions [a) 0.1 1, b) 0.08 m 3 m' 3 ] in the Southern Great Plains using the a) Noah and b) CLM 
LSMs with the YSU (red solid), and MYJ (green solid), and MRF (blue solid) PBL schemes. 
Also shown are the surface (V S f C ) and entrainment (V en t) vectors (dashed lines), surface (/fy) and 
entrainment ( (3 ent ) Bowen ratio values, and heat (Ah) and moisture (A/ e ) entrainment ratios. 

Figure 4: Same as Figure 3, but for intermediate soil moisture conditions [0.18 m 3 m' 3 ]. 

Figure 5: Same as Figure 3, but for wet soil moisture conditions [a) 0.32 and b) 0.40 m 3 m' 3 ]. 
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Figure 6: Relationship between latent ( LE ) and sensible (//) heat fluxes from the surface (■), 
entrainment (A), and their sums (0) as simulated by LIS-WRF using the six PBL-LSM 
combinations for dry, intermediate, and wet soil conditions and derived using the mixing 
diagrams in Figs. 3-5 and values given in Table 1. The dashed lines represent a theoretical 
constant surface (H + LE) net radiation of 500 Wm" 2 for comparison with the simulations. 

Figure 7: Relationship of evaporative fraction ( LE s f c /(LE s f c + H s f c )) to maximum PBL height as 
simulated by the Noah (o) and CLM (•) LSMs coupled with the YSU (red), MYJ (green), and 
MRF (blue) PBL schemes in LIS-WRF on June 12 at the dry, intermediate, and wet soil 
locations shown in Figs. 3-5. 

Figure 8: Diurnal co-evolution (12-00Z) of 2m-specific humidity ( Lq ) and 2m-potential 

temperature ( C p 6 ) on 6 June 2002 as simulated by LIS-WRF for the ARM-SGP CF at Lamont, 
OK using the a) Noah and b) CLM models and PBL combinations with the associated surface 
and entrainment vectors and derived metrics. Also overlain are observations from CF and 
metrics calculated from surface meteorology, flux, and profile measurements (black). 

Figure 9: Same as Figure 8 for the ARM-SGP Extended Facility at Plevna, KS. 

Figure 10: Same as Fig. 6, but for the a) E13 and b) E4 sites on 6 June 2002 (Figs. 8 and 9) 
along with observations (open black). The dashed line represents values of constant available 
energy (H s f c + LE s f c ) equal to that observed. 

Figure 11: Relationship of evaporative fraction (LE s f c /(LE s f c + H s f c )) to maximum PBL height as 
simulated by the Noah (o) and CLM (•) LSMs coupled with the YSU (red), MYJ (green), and 
MRF (blue) PBL schemes in LIS-WRF on June 6 at the E13 and E4 sites shown in Figs. 8 and 9. 
Figure 12: Diurnal evolution (12-00Z; solid line) of 2m-potential temperature (C p 9) vs. 2m- 
specific humidity (Lq) from a representative day during June 2002 in the U.S. Southern Great 
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Plains. The annotations on the plots depict the vector component contributions of surface and 
entrainment fluxes, and the addition of a vector due to the horizontal advection of heat and 
moisture (V a dv). 

Figure 13: Mixing diagrams generated from LIS-WRF simulations using Noah and CLM for 
wet soil moisture conditions (see Fig. 5) with the addition of the horizontal advection vector 
(V a dv) and Bowen ratio (/3 w v ). The arrows indicate the direction of the advective fluxes. 

Figure 14: Total mean fluxes of heat (H to t) and moisture ( LE tot ) and their component fluxes 
from the land surface (H s f c , LE s f c ), advection (H a dv, LE a dv ), and entrainment (H ent , LE ent ) derived 
from the a) Noah and b) CLM mixing diagrams both without (from Fig. 5) and with (Fig. 13) 
advection, and c) plotted with advection as in Fig. 6c. 
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Figure 1: a) Diurnal evolution (12-00Z; solid line) of 2m-potential temperature (C p 6) vs. 2m- 
specilic humidity (Lq) from a representative day during June 2002 in the U.S. Southern Great 
Plains as simulated by a coupled mesoscale model. The annotations on the plots depict the 
vector component contributions of surface and entrainment fluxes that are obtained using the 
mixing diagram approach, b) 12-hour time series of 2m potential temperature (0) and 2m 
specific humidity ( q ) corresponding to Fig. la. c) 12-hour time series of PBL height (PBLH) 
and surface sensible heat flux (H s f c ) which are used in calculating the vector and flux components 
in Fig. la. 
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Figure 2: Soil moisture (volumetric: m 3 m' 3 ) in the upper (0-10cm) layer valid at 12Z on 12 
June 2002 as simulated from a 2.5 year spinup of the Noah model over the 1-km LIS-WRF 
domain in the SGP for the IHOP-02 experiment. The ARM-SGP Central Facility (CF) at 
Lamont, OK is also shown along with the dry, intermediate, and wet soil locations presented in 
Figs. 3-5. 
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Figure 3: Diurnal co-evolution (12-00Z) of 2m-specific humidity (Lq) and 2m-potential 

temperature (C p 6) on 12 June 2002 as simulated by LIS-WRF during dry soil moisture 
conditions [a) 0.1 1, b) 0.08 m 3 m°] in the Southern Great Plains using the a) Noah and b) CLM 
LSMs with the YSU (red solid), and MYJ (green solid), and MRF (blue solid) PBL schemes. 
Also shown are the surface (V S f C ) and entrainment (V en t) vectors (dashed lines), surface (/3 S / C ) and 
entrainment (Pent) Bowen ratio values, and heat (A/,) and moisture (A/ e ) entrainment ratios. 
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Figure 4: Same as Figure 3, but for intermediate soil moisture conditions [0.18 m 3 
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Figure 5: Same as Figure 3, but for wet soil moisture conditions [a) 0.32 and b) 0.40 m 3 
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Figure 6: Relationship between latent ( LE ) and sensible (If) heat fluxes from the surface (■), 
entrainment (A), and their sums (0) as simulated by LIS-WRF using the six PBL-LSM 
combinations for dry, intermediate, and wet soil conditions and derived using the mixing 
diagrams in Figs. 3-5 and values given in Table 1. The dashed lines represent a theoretical 
constant surface (H + LE) net radiation of 500 Wm" 2 for comparison with the simulations. 
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Figure 7: Relationship of evaporative fraction ( LE s f c l(LE s f c + 1 1 s fc)) to maximum PBL height as 
simulated by the Noah (o) and CLM (•) LSMs coupled with the YSU (red), MYJ (green), and 
MRF (blue) PBL schemes in LIS-WRF on June 12 at the dry, intermediate, and wet soil 
locations shown in Figs. 3-5. 
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Figure 8: Diurnal co-evolution (12-00Z) of 2m-specific humidity (Lq) and 2m-potential 

temperature ( C P B ) on 6 June 2002 as simulated by LIS-WRF for the ARM-SGP CF at Lamont, 
OK using the a) Noah and b) CLM models and PBL combinations with the associated surface 
and entrainment vectors and derived metrics. Also overlain are observations from CF and 
metrics calculated from surface meteorology, flux, and profile measurements (black). 
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Figure 9: Same as Figure 8 for the ARM-SGP Extended Facility at Plevna, KS. 
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Figure 10: Same as Fig. 6, but for the a) E13 and b) E4 sites on 6 June 2002 (Figs. 8 and 9) 
along with observations (open black). The dashed line represents values of constant available 
energy (H s j c + LE s f c ) equal to that observed. 
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Figure 11: Relationship of evaporative fraction (LE s f c /(LE s f c + H s f c )) to maximum PBL height as 
simulated by the Noah (o) and CLM (•) LSMs coupled with the YSU (red), MYJ (green), and 
MRF (blue) PBL schemes in LIS-WRF on June 6 at the E13 and E4 sites shown in Figs. 8 and 9. 
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Figure 12: Diurnal evolution (12-00Z; solid line) of 2m-potential temperature (C p 6) vs. 2m- 
specific humidity (Lq) from a representative day during June 2002 in the U.S. Southern Great 
Plains. The annotations on the plots depict the vector component contributions of surface and 
entrainment fluxes, and the addition of a vector due to the horizontal advection of heat and 
moisture (V a dv)- 
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Figure 13: Mixing diagrams generated from LIS-WRF simulations using Noah and CLM for 
wet soil moisture conditions (see Fig. 5) with the addition of the horizontal advection vector 
(V adv) and Bowen ratio ((3adv)- The arrows indicate the direction of the advective fluxes. 
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Figure 14: Total mean fluxes of heat (H to i) and moisture ( LE tot ) and their component fluxes 
from the land surface {H s f c , LE s f c ), advection (H a dv, LE a dv), and entrainment (H enh LE ent ) derived 
from the a) Noah and b) CLM mixing diagrams both without (from Fig. 5) and with (Fig. 13) 
advection, and c) plotted with advection as in Fig. 6c. 


Table 1: Mean fluxes of heat (. H tot ) and moisture ( LE tot ) and their component fluxes from the 
land surface (H s f c , LE s f c ) and entrainment (H ent , LE ent ) as simulated by LIS-WRF using the Noah 
and CLM LSMs with the YSU, MYJ, and MRF PBL schemes. The flux values (in Win" 2 ) were 
derived using the mixing diagram theory and surface and entrainment flux vectors for dry, 
intermediate, and wet soil conditions depicted in Figs. 3-5. 


DRY 

H tot 

H s fc 

H ent 

LEtot 

■zm 

LEent 

Noah- YSU 

1071.7 

302.3 

769.4 

-2428.2 

96.2 

-2524.4 

Noah-MYJ 

954.2 

356.0 

598.2 

-1664.1 

94.7 

-1758.9 

Noah-MRF 

854.2 

303.6 

550.6 

-1659.1 

96.6 

-1755.7 

CLM- YSU 

1517.3 

347.7 

1169.6 

-2607.0 

41.8 

-2648.8 

CLM-MYJ 

1059.7 

346.5 

713.2 

-1381.5 

39.0 

-1420.5 

CLM-MRF 

1476.7 

339.1 

1137.6 

-2281.5 

41.3 

-2322.7 

INT 

H tot 

jggg| 

H ent 

LEtot 


LEent 

Noah- YSU 

311.7 

309.3 

2.4 

-115.2 

114.2 

-229.4 

Noah-MYJ 

374.0 

343.4 

30.6 

-119.0 

110.8 

-229.8 

Noah-MRF 

270.3 

311.3 

-40.9 

-100.1 

113.3 

-213.4 

CLM- YSU 

330.2 

281.2 

48.9 

-120.2 

34.9 

-155.1 

CLM-MYJ 

433.2 

287.9 

145.4 

-118.8 

34.1 

-152.9 

CLM-MRF 

338.0 

283.8 

54.2 

-117.9 

34.6 

-152.6 

WET 

H tot 

H s fc 

H ent 

LEtot 


LEent 

Noah-YSU 

141.8 

113.2 

28.6 

120.2 

450.8 

-330.6 

Noah-MYJ 

173.6 

133.5 

40.1 

59.5 

425.2 

-365.7 

Noah-MRF 

125.1 

115.8 

9.4 

71.4 

447.7 

-376.2 

CLM- YSU 

160.1 

-24.2 

184.3 

-45.0 

522.6 

-567.6 

CLM-MYJ 

243.2 

20.2 

222.9 

-7.9 

481.0 

-488.9 

CLM-MRF 

140.0 

-17.6 

157.6 

13.7 

516.5 

-502.8 


















































































































































